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Abstract 



We calculate the magnetic response of ensembles of small two-dimensional 
structures at finite temperatures. Using semiclassical methods and numerical 
calculation we demonstrate that only short classical trajectories are relevant. 
The magnetic susceptibility is enhanced in regular systems, where these tra- 
jectories appear in families. For ensembles of squares we obtain a large para- 
magnetic susceptibility, in good agreement with recent measurements in the 
ballistic regime. 
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A free electron gas at temperature T and magnetic field H such that ksT ^> hw (w = 
eH/ mc) exhibits a small diamagnetic response jl] . This behavior persists when the electrons 
are placed in periodic or weak-disorder potentials 0. When the system is constrained to 
a finite volume the confining energy appears as a relevant scale giving rise to finite-size 
corrections to the Landau susceptibility. These corrections have been the object of several 
theoretical studies in the last few years for the case of clean [||,|| and disordered systems, 
and received renewed interest with recent experiments of Levy et al. || : Measurements on 
an ensemble of 10 5 microscopic, phase-coherent, ballistic squares lithographically defined 
on a high mobility GaAs heterojunction yielded a large paramagnetic susceptibility at zero 
field, decreasing on the scale of approximately one flux quantum through each square. These 
experiments have been important in orienting the theoretical studies towards the physically 
relevant questions associated with the magnetic response of small systems. In particular, 
the role of finite temperature and the necessity of distinguishing individual from ensemble 
measurements appear as important ingredients that have been overlooked in some of the 
theoretical literature. 

In this letter we calculate the magnetic susceptibility of noninteracting electrons at finite 
temperatures in clean regular geometries (i.e. squares and circles) for individual systems as 
well as for ensembles. We use semiclassical approximation and classical perturbation theory 
since the magnetic fields involved are not big enough to modify the classical trajectories 
significantly. We explore the validity of our assumptions and analytical results with numer- 
ical calculations. We compare the results obtained for ensembles of regular structures with 
those of chaotic billiards, finding important quantitative differences. We show that, within a 
semiclassical approach, finite temperature induces a cut off on the classical trajectories con- 
sidered, and therefore clean systems can provide a good description of the ballistic regime. 
This is the case for the experimental conditions of Ref. || and therefore our model yields 
results in good agreement with the measurements. 

We consider an ensemble of isolated two-dimensional systems at temperature T. For 
each member of the ensemble (with N electrons and area V) the magnetic susceptibility \ 
is given by the change of the free energy F(T, N, H) under the effect of a magnetic field, 

The necessity of using the canonical ensemble for isolated mesoscopic systems, and the 
physical differences with the grand-canonical ensemble (GCE, where the system responds to 
the magnetic field with a fixed chemical potential fi), are some of the important concepts that 
recently emerged in the context of persistent currents H . On the other hand, calculations in 
the GCE are more easily performed due to the simple form of the thermodynamic potential 

n(T,fi,H) = ~J dEp(E) In (1 + exp [/?(//-£)]) , (2) 

in terms of the density of states p{E) = — (2/tt) Im g{E). The factor of 2 takes into account 
spin degeneracy, (3 = 1/ksT, and g(E) is the trace of the Green function GE(r',r)i i-e. 

g(E) = JdrG E (r,r). (3) 
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Separating p into a mean (Weyl) part, which is field independent, and an oscillating part, 
p(E) = p°(E) + p osc (E), we define a mean chemical potential p° from N = J dE p{E) f (E — 
p) = J dE p° (E) f (E — p°) . (/ is the Fermi-Dirac distribution function.) Considering that 
p osc <C p°, it has been shown || 

F(N) = F° + AF« + AF® , (4) 

where F° = p°N+n°(p°) and AF^ = fi osc (/i°). We define and fi osc by using respectively 
p° and p osc instead of p in Eq. (E3). The second-order term is [0] 



AF^ ' 



2p°(/i°) 



| dE p^F) /(£-// 



(5) 



F° is field independent and does not contribute to x- AF^ gives the susceptibility in a 
GCE with chemical potential p°. In disordered systems it vanishes under impurity average, 
and we will show that it is also the case within the averages of our semiclassical model. 

We will calculate p osc from the semiclassical expansion of the Green function. Except for 
a logarithmic singularity when r' — > r, which yields the smooth part p° of p, the semiclassical 



Green function has the generic form [10 



G-(r', r) = E A exp [1 (s t - (rj t - ±) |)] , (6) 

where the sum runs over all classical trajectories t joining r to r' at energy E. S t is the 
action integral along the trajectory. For billiards without magnetic field we simply have 
S t = hkL t , k = y2mE/h and L t is the length of the trajectory. The amplitude D t takes care 
of the classical probability conservation and r] t is the Maslov index. 

Within our semiclassical approach, the free energy corrections are given as sums over clas- 
sical trajectories, each term being the convolution in energy of the semiclassical contribution 
(oscillating as kL t ) with the Fermi factor (smooth on the scale of 0). It can be shown (Tl| 



that the T = contribution to AF^ of a trajectory is reduced by a temperature-dependent 
factor R(T) = (L t /L c ) sinh -1 (L t /L c ), with L c = h 2 k F (3 / '(irm). A factor of R 2 (T) is needed 
for AF^ . For long trajectories and high temperatures, R(T) results in an exponential 
suppression and therefore the fluctuating part of the free energy, and Xi are dominated by 
trajectories with L t < L c , which will be the only ones considered in our analysis. (We will 
not write R(T) and R 2 (T) in the equations that follow.) 

The standard route to obtain p osc from G s £ is to evaluate the integral of Eq. (|3]) by 
stationary-phase approximation. This selects the trajectories which are not only closed in 
configuration space (r' = r), but also closed in phase space (p' = p), i.e. periodic orbits. 
When these latter are [well] isolated the Gutzwiller Trace Formula [Joj is obtained. For 



integrable systems, periodic orbits come in continuous families corresponding to the rational 
invariant tori (Berry- Tabor Trace Formula [0). The difficulty in following this approach 
in our case stems from the fact that in calculating \ for small fields, one is actually looking 
at the effect of a small perturbation on rational tori. The Poincare-Birkhoff theorem states 
that, as soon as the field is turned on, generically (the circular billiard being a notable 
exception) all rational tori (i.e. all families of periodic orbits) are instantaneously broken, 
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leaving only two (one stable, one unstable) isolated periodic orbits. On the one hand, the 
physical effect which generates x is the breaking of the rational tori, so that just ignoring 
this, i.e. using the Berry- Tabor Formula, is certainly inadequate. On the other hand, for 
H— >0, the remaining orbits are not sufficiently well isolated to apply the Gutzwiller Trace 
Formula. Therefore, a uniform treatment of the perturbing field is needed, where not only 
orbits that are closed in phase space are taken into account, but also trajectories closed in 
configuration space which can be traced to periodic orbits when H—>-0. 

In squares (of side a), due to the simplicity of the geometry, such a uniform treatment is 
possible since we can perform the corresponding integrals exactly. For H = 0, rj t is twice the 
number of reflections, and D t = a/L l J 2 with a = — 7r(2m) 3//4 /[(27r/i) 3 / 2 i? 1 / 4 ]. One way to 
obtain this result is to use the method of images and express Ge in terms of the free Green 



function G° E as 



GeW, r) = G|(r>)+5>G°(r' 4l r) l (7) 



where the r- are all the mirror images of r' by any combination of symmetries across the 
sides of the square, and q = ±1 depending on the number of symmetries needed to map 
r' on r^. The long-range asymptotic behavior of the two-dimensional free Green function 
G^(r-,r) ~ aexp [-i ( A: | r ^ — r| — 7r/4)]/|r- — r] 1 / 2 can be used for the images [13 . 



For sufficiently weak magnetic fields, one may keep in Eq. (0) the zero-order approxima- 
tion for D t , and use the first-order correction 5S to the action. For a closed orbit enclosing 
an algebraic area A, classical perturbation theory yields SS = (e/c)HA for low fields and 
high energies, such that the cyclotron radius of the electrons is much larger than the typical 
size of the structure. 

We now specify the contribution p n [to p osc ] of the family of closed trajectories which, for 
H —>0, tends to the family of shortest periodic orbits with non-zero enclosed area. We note it 
(1,1) since the trajectories bounce once on each side of the square (upper inset, Fig. |IJ). Their 
length is Lxi = 2\[2a. This family gives the main contribution to the experiment of Ref. || 
since L c « 2a at T = 40mK. The contribution of other families is obtained essentially 
in the same way. However, strong flux cancellation occuring for other primitive orbits 
makes their contribution irrelevant in the case of the square, even for very low temperatures 
|lll , |14|| . Using as space coordinates x , which labels the trajectory, s the distance along the 



trajectory, and the index e = ±1 specifying the sense of motion the area is simply 
A e (x ) = e2x (a — x ). Inserting A in Eq. (|3]) we have pn(H) = p n (H=o)C(H), where 
p 11 (i/=o) = —8a 2 a sin (kLn + TT /A) / Ll{ 2 is the unperturbed contribution and 

ra /2g \ 1 

C(H) = / da; cos ( —Hxo(a — xq) ) = — [cos(7ry9)C( v /7r^) + sin(7r<p)S( v /7r(/?)] . (8) 
Jo \nc ) v2(p 



2if 

C and S are respectively the cosine and sine Fresnel integrals, and if = $/$o is the total 
flux $ = Ha 2 inside the square measured in units of $o — hc/e. For (p > 1 the Fresnel 
integrals can be replaced by their asymptotic value 1/2, which amounts to evaluating C((p) 
by stationary phase, i.e. C s (ip) = cos(7r ip + it/ '4) / This expression however diverges for 
H^O, while C(0) = 1. 

The contribution of the (1,1) family to AF^ yields, in leading-order in k F a 
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(M) 3/2 sin k F L u + -)-—. (9) 



Therefore, the susceptibility of a given square can be paramagnetic or diamagnetic (Fig. [1]) 
and its typical magnitude is much larger than Xl, with — \l = — e 2 / (12irmc 2 ) being the two- 
dimensional Landau susceptibility. Clearly, x vanishes under average if the dispersion of 
k F a across the ensemble is of the order of 2tt. The average \ is then given by the contribution 
of the (1,1) family to AF® 

M = _ 3 it o ^ . (10) 

The average susceptibility (solid line, Fig. ||) is paramagnetic at H = and for low fields it 
oscillates with an overall decay of l/ip. The divergent susceptibility obtained from C s (dotted 
line) provides a good description of x f° r f > 1- For ensembles with a wide distribution 
of lengths a (in Ref. || the dispersion in size across the array is estimated between 10 
and 30%) a second average in d 2 C 2 /d<£> 2 should be performed. (Since the scale of variation 
of C with a is much slower than that of sin (k F Ln) we can effectively separate the two 
averages.) The low-field oscillations of (x) are suppressed under the second average (done 
for a gaussian distribution with a 30% dispersion, dashed line) , while the zero-field behavior 
remains unchanged. 

We checked the semiclassical results calculating the first 1500 eigenenergies of a square 
in a magnetic field by direct diagonalization. At T = the free energy reduces to the total 
energy and x is dominated by big paramagnetic singularities at the level-crossings of states 
belonging to different symmetry classes and at small avoided-crossings between states with 
the same symmetry ||. These peaks are compensated once the next state is considered, and 
therefore disappear at finite temperature where the occupation of nearly degenerate states 
becomes almost the same. Temperature regularizes the T = singular behavior, and of 
course, describes the physical situation. We include it by calculating the partition function 



Z = exp [— (3F] from a recursive algorithm |T6| , pT|| . The results for individual squares are 



in excellent agreement with Eq. (Q), the oscillations as a function of k F Ln (and <p) clearly 
shown in Fig. |l|. The average values also agree with our analytical findings (Fig. |^). 

Ref. yielded a paramagnetic susceptibility at if = with a value of approximately 100 
(whithin a factor of 4) in units of Xl- For the two electron densities n s = 10 11 and 3x10 



ii 



cm 



of the experiment, the factor 4\/2/(5ii)k F a from Eq. (|T(]) gives respectively a susceptibility 
of 130 and 220, that when temperature is considered (through R 2 (T = AOmK)) become 60 
and 170, in good agreement with the measurements. The field scale for the decrease of 
(x(h=o)) is of the order of one flux quantum through each square, in reasonable agreement 
with our theoretical findings. The temperature scale for the decrease of the susceptibility 
was identified as given by the inverse time-of-flight Vp/a, which is the same scale L t /L c that 
we find. 

Squares constitute a generic example of an integrable system perturbed by a magnetic 
field. It is interesting to compare our results with two extreme cases: circles (which remain 
integrable under the perturbation) and completely chaotic systems. Expressing the hamil- 
tonian of a circle (of radius a) in action-angle variables [17]], p osc can be written as a sum 
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over families of periodic trajectories | 12| . Within our finite-temperature approach we restrict 



ourselves to the shortest ones, the whispering-gallery trajectories who turn only once around 
the circle in coming to the initial point after M bounces. Their contribution to p osc is 

oo / cH \ 

p wg (H) = Pm( h =°) cos — A m . (11) 

M=3 ^ nC ' 

p M ( H =o) = V8mLll 2 /( v ^h 2 k 1 / 2 M 2 ) sin (kL M + ir/A-3irM/2) and the length of the M th tra- 
jectory is L u = 2Ma sin (jr/M), while the enclosed area is A M = (Ma 2 /2) sin (2ir/M). 

oscillates as a function of k F a with an amplitude proportional to (k F a) 3 / 2 (consistent 
with Ref. and vanishes under ensemble average. {x( H=0 )) /Xl = (24/7i)k F aC , with 

C = J2m=3 s i n {k/M) sin 2 (2ir/M)/M ~ 0.7. The sums over M are rapidly convergent, indi- 
cating the dominance of the first few periodic orbits. 

Squares and circles give the same dependence on k F a for x an d (x)- This generic 
behavior for integrable systems can be traced to the k~ 1 / 2 dependence of the contribution to 
p osc most sensitive to the magnetic field. The numerical prefactors depend on the specifics 
of the geometry. The main contribution to x at if = comes from interference between pairs 
of time-reversed trajectories. In the circle all periodic orbits within a family have the same 
area, while for squares the dominant family (1,1) also includes periodic orbits with small 
enclosed area. This difference results in a larger value of (x{h=o)) (by a factor of 10) for the 
circle and the absence of the 1/tp damping of the low field oscillations. 

For chaotic systems (of typical length a) with hyperbolic periodic orbits the Gutzwiller 
Trace Formula provides the appropriate path to calculate p osc . When only a few short 
periodic orbits are important, x can have any sign and its magnitude is of the order of 
k F aXh im . Extending this analysis to the case of an ensemble of chaotic systems we obtain 
(x) Xl- The individual x are larger, by a factor [k F a) 1 ^ 2 in regular geometries than in 



chaotic systems [20]. For (x) the difference is even larger (factor k F a). This is due to the large 
oscillations of p in regular systems induced by families of periodic trajectories. The different 
magnetic response according to the geometry does not arise as a long-time property (linear 
vs. exponential trajectory divergences) but as a short-time property (family of trajectories 
vs. isolated trajectories). This assures that small variations in the geometry of the clean 
systems that we have considered will not be relevant. 

We believe that measuring the susceptibility in different geometries will be of high interest 
in order to understand the applicability of simple noninteracting semiclassical models to 
actual microstructures. 

We acknowledge helpful discussions with H. Baranger, O. Bohigas, Y. Gefen, 
M. Gutzwiller, L. Levy, N. Pavloff, R. Prange, B. Shapiro, and H. Weidenmiiller. KR 
acknowledges financial support by the Humboldt foundation. 
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FIGURES 



FIG. 1. Magnetic susceptibiltiy of a square as a function of k F a from numerical calcula- 
tions at zero field and a temperature equal to 10 level-spacings. The number of electrons is 
N = (k F a) 2 /(27r). The dashed line is the envelope of the oscillations (in k F L\i) of our semiclas- 
sical approximation with the temperature correction factor R{T). The period ir/y/2 indicates the 
dominance of the shortest periodic orbits enclosing non-zero area with length L\\ = 2y2a (upper 
inset). Lower inset: amplitude of the oscillations (in k F Ln) of x as a function of the flux through 
the sample from Eq. (g) (dashed) and numerics (solid). 

FIG. 2. Average magnetic susceptibiltiy for an ensemble of squares from Eq. ( |i~0| ) (solid) and 
from the stationary-phase integration C s (dotted). Dashed: average over an ensemble with a large 
dispersion of sizes (see text), Thick dashed: average from numerics. Inset: average susceptibility 
as a function of k F a for various temperatures (4,6 and 10 level spacings) and a flux ip = 0.15, from 
Eq. ( |i~0D (dashed) and numerics (solid). 



